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We present a comprehensive framework for treating the nonlinear interaction of few-cycle pulses using an 
envelope description that goes beyond the traditional SVEA method. This is applied to a range of simulations 
that demonstrate how the effect of a x'^' nonlinearity differs between the many-cycle and few-cycle cases. Our 
approach, which includes diffraction, dispersion, multiple fields, and a wide range of nonlinearities, builds upon 
the work of Brabec and Krausz| 1] and Porras|2]. No approximations are made until the final stage when a 
particular problem is considered. 

The original version (vl) of this arXiv paper is close to the published Phys.Rev.A. version, and much smaller 
in size. 

PACS numbers: 42.65.Re,42.65.Yj,42.65.-k,31.15.-p 



NOTE: this is a longer version of a paper published as 
Phys. Rev. A 67, 023813 (2003). It therefore contains the 
same text and graphics, but has additional remarks and many 
extra figures, most notably more complete sets of simulation 
results. This has led to some repetitions in the text, and also 
some of the extra figures have artifacts: for example, the enve- 
lope phase graphs when the envelope has a small magnitude, 
and some of the figures showing the phase drift have spikes 
where a stable value could not be found. Further, some fig- 
ures pass without much comment. 



I. INTRODUCTION 

The analysis of optical pulse propagation traditionally in- 
volves describing a pulse in terms of a complex field envelope, 
while neglecting the underlying rapid oscillations at its carrier 
frequency. The resulting "slowly varying envelope approxi- 
mation" (SVEA) (see e.g. |3|), which reduces second order 
differential equations to first order, is valid when the enve- 
lope encompasses many cycles of the optical field and varies 
slowly. The alternative approach is to solve Maxwell's equa- 
tions numerically (see e.g. |ll|4|]), which is more general but 
involves greater computational effort, and lacks the intuitive 
picture of a pulse "envelope". 

For example, optical parametric oscillators (OPOs) based 
on aperiodically-poled lithium niobate (APPLN) have gener- 
ated 53 fs idler pulses at 3/jm that are nearly transform limited, 
and contain only five optical cycles |5]; laser pulses with less 
than three optical cycles have been generated in other contexts 
10]. Under these circumstances, the validity of the slowly- 
varying envelope approximation is clearly open to question. 

Brabec and Krausz |1] derived corrections to the SVEA, 
which they included in their "slowly evolving wave approx- 
imation" (SEWA). This enabled the few-cycle regime to be 
modelled with improved accuracy, and the SEWA has subse- 
quently been applied in different situations, including ultra- 
short IR laser pulses in fused silica the filamentation of 
ultra-short laser pulses in air |9|, and even in microstructured 
optical fibres lUOIl . Later, Porras |2] proposed a slightly dif- 
ferent "slowly evolving envelope approximation" (SEEA) that 



included corrections for the transverse behavior of the field. 

Here we use a field envelope approach to simulate the prop- 
agation of ultrashort pulses in a x'^'medium. The novelty is 
that we (a) derive a more general form than that of Brabec 
and Krausz, called the "generalised few-cycle envelope ap- 
proximation"(GFEA); and (b) apply it to both optical (non- 
degenerate) parametric amplification (NPA), and the optical 
parametric oscillator (OPO). A shorter version of this paper 
has been published in Phys. Rev. A [ll]. The only previous 
multiple field app lication of this kind of result was for four 
wave mixing lUlll . 

We compare the SEWA/SEEA equations to our own (sec- 
tion inj, and explain the differences and subsequent adjust- 
ments to the necessary approximations. This theory enables 
us to rigorously study what combination of approximations af- 
fords the most efficient method for treating a given nonlinear 
interaction involving few-cycle pulses. Next (section HIH we 
discuss the x'^' nonlinearity and a scaling scheme designed 
to reveal the few-cycle effects. Then we compare the SVEA 
predictions to the few cycle GFEA theory using idealised sit- 
uations (section nVl and more realistic OPO models (section 
|yj. Finally, we present our conclusions (section lvH . 



II. THEORY 

This section contains a summary of a complete rederivation 
llisll of a Brabec & Krausz style theory, which yields an evolu- 
tion equation for an envelope description of pulse propagation 
in the few-cycle regime. Our result is more complicated than 
the SEWA equation but reduces to it in the appropriate 
limits; it also explains the slight differences between their re- 
sult and that of Porras 01 . 

Following Brabec-Krausz, we consider the case of small 
transverse inhomogeneities of the polarization, and so start 
with the three dimensional wave equation 

{dl + V\)E{rf)^^d^ f dt'e{t~t')E{r,t') 

C J — oo 



47t 



d}P„i{rf). 



(1) 
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Here is the transverse Laplace operator, 3a is shorthand 
notation for d/da, e(f) = (271)"' /r„t/a)e(a))e'™, e(co) = 1 + 
47tx(co), and x{(o) is the linear electric susceptibility. The 
electric field E propagates along the z direction. Both E and 
the nonlinear polarization P„i are polarized parallel to the x 
axis. 

We can transform eqn. Q into frequency space in order to 
expand e(co) in powers of CO, thus enabling us to treat the ma- 
terial parameters as a power series which we can truncate to 
an approriate order However for simplicity it is better to ex- 
pand k about a suitable coq instead. Using e(co) — c^A;(a))^/co^, 
it follows that 



Y„(a)-cooj 



y„ 3m^(co)|„g = p„ + ia„;p„,a„ e : 



(2) 



We can now write the frequency space version of eqn. ([0 



and eqn. (|5j now becomes 



(ipo+M5-cooPiax)'+vi 



n=0 



xA(r^,^,T) + ^ (1 + id,fB{r^,^,T,A) = 0. (8) 



For convenience we also introduce the dimensionless 
ratio of phase and group velocities a = COopi/po = 
(coo/po) /(I /Pi) = Vf/vg, and use the fact that the refractive 
index at coo is no = cpo/coo. We also define a dispersion term 
D in a similar way to Brabec-Krausz, but instead use a scaled 
(dimensionless) version D' ~ (coo/Po)^ in following equa- 
tions so that 



b' 



COo 

Po 



E^^K)" 



n=l 



(9) 



(a2 + Vi)£(r,f)- 
= —djPni{r,t). 
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("Y„(ar + iCOo)" 



E{r,t) 



(3) 



We introduce an envelope and carrier form for the field in 
the usual way, using r = {r^,z), so that 



E[r,t) = A(rx,Z,f)e'^ +A*(rx,z,f )e" 



(4) 



and similarly P„i{r,t) = B{r^,z,t;A)e''^ + B*{r±,z,t;A)e^''^. 
The symbol E — poz — coo? + V|/o is introduced as a convenient 
shorthand for the argument of the exponential. With these 
envelope-carrier substitutions, the equation of motion can be 
written as 



([iPo + 3/ + vi)A(?x,z,r) 

to «' V«o J \ 

^ 1 + — a, B{r^,z,t;A). 

\ f\\r\ I 



(5) 



Eqn. (|5} has the opposite sign on the RHS to Brabec & 
Krausz's eqn. (2), but agreement is recovered later in eqn. 
(O- 

As is usual, we introduce co-moving variables 



T = COo (f - Piz) , 3r = «o3x, 
^ = Poz, 3, = Po35-cooPi3t, 



(6) 
(7) 



Hence we get 



ii 

ia{id^) - — + iD' 
Po 



la (id-c) - — + iD' 
po 



2% 

+ —il+id,fBir^,^,V,A). (10) 
'«() 

This form can be rearranged without approximation to 



+ —i^. ^B{r^,^,z;A) 



Tr 



nl (l + !a8x) 



1 + lOdr 



(11) 



where 



Tr 



iq^ 2 , ' A «0 



lb' 



A(rx,^,T). (12) 



Eqn. (II 1> exflcf - it contains no more approximations 
than our starting point eqn. ([0 except for the expansion of e 
in powers of CO. We recover the full field E from eqn. @ by 
recombining A and knowing the carrier. The partial deriva- 
tives (;8x) in the denominators can, if necessary, be treated by 
Fourier transforming into the conjugate frequency space (i2). 
Note that like i, £1 is scaled relative to the carrier frequency. 
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If we set Tr = 0, this gives us a generaHsed few cycle enve- 
lope (GFEA) equation, which contains the SVEA SEWA 
d], and SEEA 12] within it as special cases. Of course we 
cannot just set the Tr term to zero without some justification, 
and this is discussed below. 

The (im/nQ^ polarization term from eqn. has 
prefactors which depend on the time-derivative of the polar- 
ization, and these new terms are what add the effect of finite 
pulse lengths to the pulse evolution. Note that we can write 
this polarization term in different forms: 



5SC 



(l+ia8x) 
= (l+io3x) 



l + (l-o) 



2(3, + (1+ a) a? 



i(l-o)3x 



(13) 



With O = 1, these reduce to the 3C = 1 + '^x SEWA [1] 
form. Similarly, to first order in (a — 1 ), one can get the ^ = 
1 + ;a3x form analogous to the SEEA corrections although 
a strict expansion to first order in 3x fTsIl gives a prefactor of 
[1 + i (2 - a)3x]. Finally, for a SVEA theory, 2C = 1, since the 
3x terms are assumed to be negligible. 

The Tr term is negligible if the following conditions hold: 

Dispersion: terms in 3x can be neglected if 



Pom 

Evolution: terms in 3^ can be neglected if 

|8^A(r^,^,£2)|«|A(r^,^,n)|, 
and eqn. dlSI l only holds if, in addition. 
Diffraction: terms in can be neglected if 

(l+On)p2H;2»l, 

Nonlinearity: is "weak" if 

\B{r^.X,^;A)\ 



< U(rx,^,n) (14) 



(15) 



nl (l+oH) 
2n {l+af 



> 



(16) 



(17) 



We use Q. instead of id-^ for these conditions because con- 
ditions on the frequency components of the various terms are 
better defined than those for time derivatives. 

These conditions are the same as those required for the 
SEWA and SEEA theories, with the SVEA conditions being 
a special case gained by setting |i2| ^ 1 for the diffraction 
and nonlinearity conditions - implying that modulations in 
the envelope are so slow compared to the carrier frequency 



that they can be neglected. Note that backwardly propagating 
behaviour has not been explicitly excluded, but since it would 
appear as a modulation on the envelope A, it would be approx- 
imated away as part of the evolution condition (eqn. il5i ). 

Note that the exact solution of eqn. il 1> makes no reference 
to a particular choice of carrier phase \\to- This implies that 
once a solution for the propagation of a particular envelope 
has been obtained, the problem has in fact been solved for a 
set of pulses (and initial conditions) based on different carrier 
phases - where that set is determined by the initial envelope 
and some arbitrary choices of carrier phase G [0,27t). The 
final state is then given by the chosen combined with the 
final form of the envelope. 



III. THE x*^' NONLINEAR SYSTEM 

When modelling x'^' nonlinear systems we split the optical 
field into two or three parts, depending on whether a degener- 
ate or non-degenerate system is being treated. For example, a 
parametric amplifier would have pump, signal, and idler field 
components. We then define an envelope Aa and carrier e'"", 
~ Pa,oz — K)a,of + Wa.Q for each field component, and use a 
separate propagation equation for each. The total field is then 
the sum of these different components: 



+A,-e'^ 



' +A:e- 



(18) 



Because the wave equation eqn. Q is linear in the electric 
field, we can use eqn. ( I18t in the theory of sectionHlJ and get a 
GFEA-like equation which is rather like the sum of three sep- 
arate copies of eqn. il 1> . but with a single polarization term. 
In a x'^' medium this polarization term is proportional to the 
square of the field, and inspection of eqn. ( I18t shows that this 
can be expanded into a sum of many terms. We can then split 
the multiple-field GFEA equation into three separate GFEA- 
like equations, one for each envelope. In doing this we need 
to make sure to assign the pieces of the polarization term ap- 
propriately to a suitable envelope equation, whilst discarding 
those that are not resonant with a field carrier 

Our chosen nonlinear crystal is congruent LiNbOs, for 
which we calculate refractive index and dispersion data from 
the Sellmeier equations of Jundt 1 14]. We model the nonlin- 
ear polarization using the square of the total electric field, re- 
taining the parts resonant with our field carriers in the normal 
way. Our OPO simulations (see section |VJ assumed a pump 
frequency of 357.1 THz, with nominal signal and idler carrier 
frequencies of 257.5 THz and 99.6 THz respectively (wave- 
lengths OMOOO/jm, 1.16500A(m, 3.0110A(m). This means the 
pump pulse will travel through LiNbOs more slowly than the 
signal and it therefore needs to be injected into the crystal 
ahead of it (see fig. 0. When the two overlap, an idler is gen- 
erated by the nonlinear interaction, and the three pulses then 
continue to interact with each other as they propagate through 
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FIG. 1: Pump timing offset (see section HIH . Tiie pump pulse is 
injected into tlie crystal just before the signal pulse is reflected off 
the input mirror. The faster moving signal pulse then catches the 
pump pulse up about halfway through the crystal, and an idler pulse 
is generated. 

the crystal. Note that our ideal non-degenerate parametric am- 
plifier simulations (see section llVt use the same field frequen- 
cies, but idealise the crystal parameters by setting the group 
velocities and dispersions to zero. 



A. System Scalings 

In a typical experiment, the crystal length would be fixed, 
as would any properties defined by its design, such as peri- 
odic poling. If we were to investigate this case for a range 
of pulse durations, there would naturally be differences be- 
tween the results, even within the S VEA. For example, the rel- 
ative pulse broadening caused by travelling through a 1000/jm 
crystal is greater for a 12fs pump pulse than for a 48fs one. 
Similarly, a fixed timing offset for injection would have dif- 
ferent effects; and a fixed pump pulse power would generate 
different strengths of nonlinear interaction for different pulse 
lengths. AH these effects would confuse any attempt at a sys- 
tematic comparison of the few cycle effects in the models we 
consider. 

Therefore, in order to isolate specific few-cycle effects, we 
must scale the pump pulse FWHM 'Tp, crystal length Lq, 
pump pulse energy fU , and pump timing offset Ja in such a 
way as to ensure the effects of group velocity, pump timing, 
and nonlinearity occur in the same proportions to one another 
over the range of pulse lengths. 

We can work out an appropriate scaling by examining a 
simple version of the propagation equation (eanll 1>. where 
we write the group velocity prefactors as B i , the second order 
dispersion prefactors as B2, and the polarization terms as CA^. 
To assist us with the scaling process we also write ^ ~ r^^^', 
1 = r^^i' , and A ~ r'^A', where r is the scaling factor. Our 
simple propagation equation is therefore 

r''+fd^,A' = 1 a,, A' + r''+2«B2a?/A' + r^''CA'^ .(19) 

We can easily match the LHS term with the polarization 
term by setting f — h; but then we must choose either h = 
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g to match group velocities, or / = 2g to match the second 
order dispersion - we cannot match both. For our chosen OPO 
situation (see section|V), it is best to match the group velocity 
terms, which control how long the pump and signal pulses 
overlap - in general, the dispersion has a much smaller effect. 

We take our reference situation to be a 20nJ 24fs FWHM 
pump pulse propagating through a 500/jm crystal, with a 
pump timing offset of 48fs. For the chosen parameter scal- 
ing (f = 8 = h) 

JIp_ ^ Lc ^ ^ 20nJ 
24fs 500A(m 48fs iV ' 

This removes the gross effects caused by reducing the pulse 
duration, so that if there are no few cycle effects, and second 
order dispersion is negligible, each pulse length should give 
identical results. 

First, the crystal length Lc is scaled in direct proportion 
to the pulse length Tp. This means that the relative group 
velocities are scaled so that the pump and signal pulses spend 
the same relative time overlapping, leading to an equivalent 
nonlinear effect. 

Second, and related to the previous point, the pump timing 
offset fA is reduced in direct proportion to the crystal length 
(and hence pulse length). This means the pump-signal pulse 
cross-over occurs at the same relative point, so the effect of 
the pump-signal overlap is similar. If a 48fs pump pulse in 
a 1000/jm crystal aiTives 96ns before the signal, the crosover 
occurs at about the 500/jm half-way mark - just as for a 12fs 
pulse in a 250/jm crystal, where a f a = 24/s gives a crossover 
at the 125/jm half-way mark. 

Third, a shorter crystal reduces the effect of the nonlinear- 
ity, so the pump pulse energy is scaled up as the inverse of the 
reduction in crystal length - a half length crystal needs pulse 
energies W twice as large - a lOnJ 48fs pulse in its 1000/jm 
crystal with see the same nonlinear effect as a 20nJ 24fs pulse 
in its 500/jm crystal. Note that the pulse energy depends on 
the area of the pulse, and so depends on its length, so that 
W oc \A\. 

We could choose to make the scaling perfect, by also scal- 
ing the crystal parameters. If we scale the crystal disper- 
sion with B2 ~ the relative amount of pulse spread- 
ing changes to become the same for each simulation - e.g. if 
the 48fs pulse widens by 10% in a 1000/jm crystal, the 12fs 
will also widen by 10% in its 250/jm crystal. We did a set of 
S VEA simulations on this basis, and as expected saw identical 
pulse profiles regardless of the chosen pulse length. However, 
we chose not to use this perfect scheme for the bulk of our 
OPO simulations because it is far from being experimentally 
practical. 



IV. IDEAL PARAMETRIC INTERACTIONS 

A parametric amplifier is a single-pass device: pump and 
signal pulses are injected into one end of the nonlinear crys- 
tal, they interact within it, then exit at the far end. However, 
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FIG. 2: SHG: Second Harmonic Generation (section llV At . Tiie 
thickness of the arrows is intended to give an indication of how the 
energy of the field components changes during propagation through 
the crystal. 



because real nonlinear crystals (such as LiNbOs) tend to have 
significant dispersion, very short pulses quickly spread out, 
making them difficult to create, and reducing the few-cycle 
effects we aim to study. 

In order to demonstrate clearly the nature of few-cycle ef- 
fects in x'^' materials, in this section we investigate an ideal 
case by setting the dispersion to zero, and make the group 
velocity the same at all frequencies. This means that O = 1, 
so the "few-cycle" polarization prefactor 3(; is identical for 
both the SEWA and GFEA theories. Note that it is difficult 
to do no-dispersion simulations over long times because pulse 
self-steepening causes both the numerical integration and the 
theoretical approximations to break down. 

We inject Gaussian pump and signal pulses at exactly the 
same time (i.e. Ja ~ 0), with the same width. They then travel 
down the crystal with maximum overlap, interacting all the 
way. Other parameters are fixed by the scaling rules from 
section nil Al Further, when graphing results for the figures, 
we scale the times for each pulse length to the 6fs case (e.g. 
for a 24fs pulse, "x — 10" corresponds to 40fs), and scale the 
pulse intensities in proportion to their initial intensities. This 
means that graphs of the initial conditions for a range of pulse 
lengths would be identical. 

Finally, note that in these ideal results, the nonlinear inter- 
action is "strong", with significant transfer of energy between 
the fields. 



A. Second Harmonic Generation (SHG) 

In second harmonic generation the nonlinearity causes a 
field at one frequency (the fundamental) to generate one at 
twice that frequency (the second harmonic). In our 48fs sim- 
ulations, a lOnJ pulse is injected in at the fundamental fre- 
quency (see fig. 13. 

In Figs. I^lsjwe can see comparisons of the pulses resulting 
from a single pass through the crystal. Although the plotted 
(and scaled) \A\^ of the pulses are rather similar; the GFEA 
theory gives noticeably different results for the shortest pulses, 
unlike the SVEA theory. Note that the vertical scales on the 
two figures (not shown) are rather different. 

However, this apparent similarity in |A|^ is a little mislead- 
ing because the underlying phase profiles of the envelopes are 
very different, as can be seen in (see fig. |6j. More features, 
primarily extra oscillations, do develop at shorter pulse dura- 
tions, but generally the comparisons remain similar to those 
shown. 




FIG. 3: SHG: Scaled pulse envelopes lA]" from the SVEA, on exit 
from the ideal dispersionless crystal. Top is fundamental, bottom is 
second harmonic. 




FIG. 4: SHG: Scaled fundamental pulse envelopes \A\ from the 
GFEA, on exit from the ideal dispersionless crystal. The SVEA re- 
sults for all pulse lengths are essentially identical to the 96fs result. 



B. Degenerate Parametric Amplification (DPA) 

A degenerate parametric amplifier uses energy from a pump 
pulse to amplify that of a signal pulse at half the frequency 
((Op = 2(x)s, see fig. Q. For our 24fs reference simulations, 
a 20nJ pulse is injected at the pump frequency, and a signal 
pulse of lOpJ is amplified; the nonlinear interaction is very 
strong. 
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FIG. 5: SHG: Scaled second harmonic pulse envelopes from the 
GFEA, see fig. |4| The S VEA results for all pulse lengths are essen- 
tially identical to the 96fs result. 
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FIG. 6: SHG: envelope-phase profiles on a scaled x, for an 18fs 

pulse duration. SVEA ( ), GFEA ( ). NB: the label "time t" 

should read "scaled t". 




FIG. 8: DPA: Scaled pulse envelopes |A| from the SVEA, on exit 
from the ideal dispersionless crystal. Top is signal, bottom is pump. 




FIG. 9: DPA: Scaled signal pulse envelopes \A\ from the GFEA, on 
exit from the ideal dispersionless crystal. The SVEA results for all 
pulse lengths are essentially identical to the 96fs result. 




FIG. 7: DPA: Degenerate Parametric Amplification (section llVBt . 
The thickness of the arrows is intended to give an indication of 
how the energy of the field components changes during propagation 
through the crystal. 



In Figs. |5J^|we can see comparisons of the pulses result- 
ing from a single pass through the crystal. Note that the inten- 
sity profiles of the pulses are very similar, but the GFEA the- 
ory gives different results to the SVEA for the shortest pulses. 

However, this apparent similarity in |Ap is a little mislead- 
ing because the underlying phase profiles of the envelopes are 
very different, as can be seen in fig. The change in phase 
of the pump pulse in both the SVEA and GFEA case is due 
to the presence of nodes that develop to either side of its cen- 
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FIG. 10: DPA: Scaled pump pulse envelopes from the GFEA, see 
fig. |9| The S VEA results for all pulse lengths are essentially identical 
to the 96fs result. 
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FIG. 11: DPA: envelope-phase profiles on a scaled X, for an 18fs 
pulse duration. SVEA ( ), GFEA ( ). . 



tre. More features, primarily extra oscillations, do develop at 
shorter pulse durations, but generally the comparisons remain 
similar. 



FIG. 12: DPD: Degenerate Parametric De-amplification (section 
IIV Ct . The thickness of the arrows is intended to give an indication of 
how the energy of the field components changes during propagation 
through the crystal. 




FIG. 13: DPD: Output SVEA signal intensities for = 3l/4 and 
= for a range of pulse durations: peak value |Ap ~ 4 x 10^; 



cause the de-amplification is very sensitive to the phase rela- 
tionship between the pulses, and the few-cycles terms modify 
the phases of the pulses. Our SVEA model of the situation 
behaves in the expected way, with the signal pulse decaying 
away to zero as it propagates through the crystal and transfers 
its energy to the pump. 

However the new terms in the GFEA propagation equa- 
tion adjust the phases of both the pump and signal fields as 
they propagate, moving them away from an exact match to 
the de-amplification criteria. As the mismatch increases, the 
de-amplification process is reversed, and instead the usual am- 




C. Degenerate Parametric De-amplification (DPD) 

This is the same situation as for degenerate parametric am- 
plification, but we have changed the relative phase between 
the pump pulse and signal so that the input signal pulse is 
de-amplified (see fig. I12> . We are interested in this case be- 



FTG. 14: DPD: Output GFEA signal intensities for = 7l/4 and 
= for a range of pulse durations, peak value 4 x 10^. The 
intensities and times are scaled in our usual way. The SVEA results 
are nearly identical for all pulse lengths, and would not show up on 
the scale of this GFEA graph: their output signal intensity consists 
of two peaks (at t ~ ±10), with a height less than l/30th of that of 
the smallest (48fs) GFEA peaks. 
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FIG. 15: DPD: Output SVEA pulse energies (Esignal = 

/ K^signai dx, arbitrary units) for a range of initial signal (envelope) 
phases (j)., and pulse lengths Tp. 



plification takes over. Fig. [H] shows the output signal pulse 
profiles on exit from the crystal, note the enormous differ- 
ence in the few-cycle dependence of the SVEA and GFEA 
theories - both in character and amplitude. The SVEA pro- 
file is just the residual input signal pulse which has not yet 
been fully "de-amplified"; this residue can also be seen as a 
shoulder on the wings of the GFEA profiles if plotted on a log 
scale. The large and very visible double peak in the GFEA 
graph is a result of the finite pulse length on the nonlinear 
polarizations term - as can be seen in eqn. jl3L this alters 
the phase profile of the nonlinear effect, and hence the pulses, 
which means the signal no longer satisfies the exact criteria 
for de-amplification. For our chosen input pulse powers, this 
effect persists for initial phase mismatches up to ^ 0.01 ra- 
dians, after which the amplification of the signal pulse has 
largely swamped the anisotropy induced by the few-cycle ef- 
fects. Of course the effect is most visible when the interaction 
has been strong enough for the input component of the signal 
pulse to (almost) completely disappear. 

In fig. [^we see how the behaviour changes both with 
pulse length and initial phase match - we see perfect de- 
amplification in the SVEA model, but the steeply sloping 
'I'signai = 7t/4 line clearly demonstrates the effects of a fi- 
nite pulse length. Of course as the (|)signai moves further 
from 7i/4, the SVEA model no longer undergoes perfect de- 
amplification, and so the difference between the two models 
becomes small. 



D. Non-degenerate Parametric Amplification (NPA) 

We consider first a non-degenerate parametric amplifier 
with pump, signal, and idler frequencies such that co,, 
G)s + ft),, and cOi ^ co,. In the 24fs reference case, the initial 
pump energy is 20nJ and the initial signal energy is lOpJ, with 
a negligible (but finite) idler. For other pulse durations, the 
energies were scaled according to eqn. ilQl . Fig. 1 171 shows 
how, according to the GFEA, the idler pulse intensity profiles 
|A, p generated in a single pass of the crystal vary with pulse 



FIG. 16: DPD: Output GFEA pulse energies (Ejignal = 

/ K^signai dx, arbitrary units) for a range of initial signal (envelope) 
phases and pulse lengths Tp. The intensities and times are scaled 
in our usual way. The SVEA results are nearly identical for all pulse 
lengths, and differ from the 48fs (i.e. logiQ(48) = 1.68) results in 
that the near (j)signal = it/4 give significantly lower energies - down 
to lO'* rather than 10^ for Tz/4. The graphed signal phases in mul- 
tiples of 3t are are 0.250, 0.251, 0.252, 0.253, 0.254, 0.256, 0.258, 
0.260, 0.275, 0.300, 0.400, 2.000 



Pump 

Signal ■ 
Idler 



FIG. 17: NPA: Non-degenerate Parametric Amplification (section 
IIV D> . The thickness of the arrows is intended to give an indication of 
how the energy of the field components changes during propagation 
through the crystal. 



duration. The profiles show little variation with pulse dura- 
tion except for the shortest pulses (x < 20), where distortion 
is evident; the signal and pump profiles show deviations of a 
comparable magnitude. 

In fig. ^]|20||^we can see equivalent comparisons for the 
signal and pump intensity profiles. 

More dramatic effects appear in the phase profiles: in fig. 
I22I the phases of the pulse envelopes at pulse durations of 
18fs and 96fs are shown with the phase distortions due to 
the finite pulse lengths (see eqn. jl3». As the pulse dura- 
tion shortens, the principal effect is to increase the magni- 
tude of the phase distortion, leaving the shape of each profile 
largely unchanged; however more complex phase oscillations 
develop for the shortest pulses. At 96fs, the profiles show 
a smaller distortion, and are tending towards the long-pulse 
SVEA limit. In this limit, the profiles are essentially flat, al- 
though the pump field develops nodes which give rise to a 
step-hke change in the phase. 



E. Non-degenerate Parametric De-ampliflcation (NPD) 

As a variant on the case just treated, signal and idler pulses 
with equal numbers of photons were injected, and the rela- 
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FIG. 18: NPA: Scaled SVEA pulse envelopes |A| on exit from the 
ideal dispersionless crystal. Top to bottom is idler, signal, pump. 



live phases of the pulses set to ensure that the signal and idler 
experience initial de-amplification (see fig. \23i . Since the 
subsequent evolution is sensitive to phase changes, and the 
finite pulse length terms in the GFEA affect the phases, this 
is an interesting situation to examine. In the SVEA, the sig- 
nal and idler decay away towards zero as the pulses propa- 
gate, so the SVEA output signal is just some residual part of 
the input. The GFEA evolution is different, as can be seen 
from eqn. ( I13> - the finite pulse lengths alter the phase pro- 
file of the nonlinearity, and hence change the evolution of the 
pulses. During an initial period of de-amplification, the pulses 
undergo a gradual phase distortion. Then, as the discrepancy 
increases, amplification takes over. In a comparison of SVEA 
and GFEA models, the effect caused by the phase distortion 




FIG. 19: NPA: Scaled GFEA idler pulse envelopes \A\ on exit 
from the ideal dispersionless crystal. The SVEA results for all pulse 
lengths are essentially identical to the 96fs result. 




FIG. 20: NPA: Scaled GFEA signal pulse envelopes, see fig. 1191 
The SVEA results for all pulse lengths are essentially identical to the 
96fs result. 



is more visible when the interaction is been strong enough 
for the input component of the signal pulse to be strongly de- 
pleted, and also is much stronger for shorter pulses 

The GFEA signal pulse profiles on exit from the crystal as 
a function of pulse duration are presented in fig. Note 
that the SVEA prediction corresponds to the long-pulse limit 
of the GFEA figure, but those limiting features are too small 
to be seen. 

The double peak is a result of the finite pulse length on 
the nonlinear polarization - the "phase twist" caused by the 
few cycle terms is odd, causing a node at T = as the field 
is amplified - at exactly T = 0, the field continues to be de- 
amplified. Note that at 18fs, for example, the output signal 
pulse amplitude is larger than the input value. 

The SVEA profile in fig. [S^is just the residual input signal 
pulse which has not yet been fully "de-amplified"; this residue 
can also be seen as a shoulder on the wings of the GFEA pro- 
files if plotted on a log scale. Fig. |25]shows the output GFEA 
signal pulse profiles on exit from the crystal, note the enor- 
mous difference in the few-cycle dependence of the SVEA 
and GFEA theories - both in character and amplitude. 
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FIG. 21: NPA: Scaled GFEA pump pulse envelopes, see fig. 1191 
The SVEA results for all pulse lengths are essentially identical to the 
96fs result. 
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FIG. 22: NPA: envelope-phase profiles for 18 and 96fs pulse dura- 
tions. Top to bottom: idler, signal, pump; SVEA ( ), GFEA 96fs 

( ), GFEA 18fs ( ). 
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FIG. 23: NPD: Non-degenerate Parametric De-amplification (sec- 
tion |l^£j. The thickness of the arrows is intended to give an in- 
dication of how the energy of the field components changes during 
propagation through the crystal. 




FIG. 24: NPD: Output SVEA signal intensities for (j)j = 3l/4 and 
= for a range of pulse durations: peak value |Ap ~ 4 x 10^). 




FIG. 25: NPD: Output GFEA signal intensities for (^s = it/4 and 
=0 for a range of pulse durations, peak value 6 x 10^. Equiv- 
alent SVEA results are very different: they are the same for all pulse 
lengths, are too small to show up on the scale of this graph (being 
~3% of the height of the 48fs GFEA peaks), and the two peaks are 
located further from the origin (at T ~ ±10). 



In fig. |5S]we see how the SVEA behaviour changes both 
with pulse length and initial phase - here we see perfect de- 
amplification. 

The GFEA output pulse energies are displayed in fig. ^\ 
which shows how the behaviour changes both with pulse 
length and initial phase. The data for (|)signai = 7t/2 demon- 
strates the effects of exact initial conditions and finite pulse 
length; maximum de-amplification occurs in the long-pulse 
(SVEA) limit. If we instead start with a signal phase slightly 
different from k/2, e.g. 0.5 l7t, the de-amplification is less ef- 
ficient and will eventually be overtaken by the amplification, 
even for the SVEA model. Consequently, comparisons for im- 
perfect initial phases are dependent on the length of the crys- 
tal. However, since we use a scaling procedure, the results still 
behave in a systematic way, even if they are not completely 
generic. 

Of course, changing other initial conditions can also disturb 
the de-amplification: e.g. different numbers of signal and idler 
photons. Although both signal and idler will initially be de- 
amplified, as they approach zero photon number, one field will 
"overshoot" the zero and be inverted. This alters the phase 
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FIG. 26: NPD: Output SVEA pulse energies (Esignal = 
/ K^signai dx, arbitrary units) for a range of initial signal (envelope) 
phases (j)., and pulse lengths Tpuise- 




FIG. 27: NPD: Output GFEA pulse energies (Esignal = 

/ K^signai dx, arbitrary units) for a range of initial signal (envelope) 
phases (|)s and pulse lengths Tpuise- The intensities and times are 
scaled in our usual way. The SVEA results are nearly identical for all 
pulse lengths, and differ from the 48fs (i.e. logig(48) = 1.68) results 
in that the near (|)signal = 7t/2 give significantly lower energies - down 
to 10"* rather than 3 x 10^ for n/2. 



relationships, and so again amplification takes over. As an 
example, simulations based on our 18fs pulses suggested that 
photon number mismatches of about one percent would not 
noticeably disrupt the appearance of either fig. ESIorET) 

These comparisons testing the sensitivity of these results to 
mismatches in the initial conditions are of course dependent 
on our chosen system parameters. For example, either more 
intense pulses or longer crystals would diminish the difference 
between SVEA and GFEA predictions; conversely, weaker 
pulses or shorter crystals would enhance them. 



V. OPTICAL PARAMETRIC OSCILLATION (OPO) 

We move on from optical parametric amplification to a 
synchronously-pumped Optical Parametric Oscillator (OPO). 
As shown in fig. |28l we considered the case of a LiNbOs 



crystal in an optical cavity with mirrors that reflect the sig- 
nal wavelength only. The oscillator is driven by a train of 
gaussian pump pulses whose periodicity closely matches the 
natural period of the cavity, and which amplify and then sus- 
tain the signal pulse confined within it. The cavity length can 
be "tuned" about exact synchronisation. The idler pulse, gen- 
erated when the signal pulse interacts with each new pump 
pulse, is transmitted through the output mirror with the pump, 
while the signal is strongly reflected. For a given set of param- 
eter values, we modelled the development of the signal pulse 
over many cavity transits until it reached a steady state. Typ- 
ically, we found that the signal stabilised in several hundred 
transits although, in a few cases, no equilibrium was achieved 
and the system oscillated indefinitely. We studied the evo- 
lution for a range of phase-matching conditions (A^ = to 
24 X 10^^//jm, and cavity length tunings (— 16fs to +16fs). 
Except where stated, the figures contain data for the perfectly 
phase matched and sychronised case. 

The steady state was determined by propagating the pulses 
through many passes of the oscillator until the modulus- 
squared of the pulse envelopes had stabilized. 



A. Dispersion Scaled OPO 

As discussed in section llll Al it is possible to make the 
SVEA propagation equation completely scale invariant if we 
modify the crystal dispersion parameters as well as the pulse 
length, energy, and crystal length. Although this is not an ex- 
perimentally achievable goal, it is instructive to look at what 
happens to the OPO output in this instance. To make auto- 
matic processing of the results more convenient, we adjusted 
the crystal lengths to be exact powers of two, which differed 
slightly from those derived from the standard reference length: 
i.e. the 24fs pulse case related to a crystal length of 512/jm, 
not 500/jm. As expected from the scaling, and shown on figs. 
EH |30| the scaled SVEA results were independent of pulse 
length. 

In contrast the GFEA pulse profiles have a more complex 
behaviour, as shown on figs. |^|32] For the chosen parame- 
ters, they have a profile which becomes increasingly distorted 
at shorter pulses; but there is also an abrupt transition to what 
is an SVEA-like pulse profile at between 36 and 48fs. 

On fig.|33] we can see this transition and the expected grad- 
ually convergence of GFEA and SVEA; the maximum differ- 
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FIG. 29: Dispersion Scaled OPO: Time domain representation of 
the S VEA modulus-squared of tiie pulse envelopes, for a range of in- 
jected pump pulse durations: 6-192fs. (bottom to top) pump, signal, 
and idler. 



X 




FIG. 30: Dispersion Scaled OPO: Time domain representation of the 
SVEA amplitude phase, for a range of injected pump pulse durations: 
6-192fs. (bottom to top) pump, signal, and idler. 



B. Scaled OPO 



ences plotted are about 10% of the pulse heights, this message 
is repeated on fig. |34] If the convergence of GFEA to SVEA 
seems slower than inutitively expected, this could be partly 
explained by noting that the effect of the extra 8x term caused 
by the non slowly-varying nature of the pulses has less dis- 
tance over which to accumulate, so a similar discrepancy be- 
tween SVEA and SEWA at e.g. 24fs and 12fs means that the 
12fs case accumulates additional evolution at twice the rate as 
at 24fs - as indeed is to be expected, since broadly speaking, 
the time derivative of a pulse envelope will double as its pulse 
length is halved. 

Note: the slight disruption to the wings of the 92fs results 
are because the time width allowed for by the simulation is 
smaller when compared to the pulse width for this particular 
case. 



The complex nature of the dynamics, which arises from 
repetitive cycling of the signal pulse in the presence of many 
interacting processes makes the isolation and analysis of few- 
cycle effects within the different models quite complicated. 
Fig. 1^3] shows intensity profiles for the pump, signal, and 
idler (bottom to top in each frame) for the SVEA (dashed 
line) and GFEA (solid line) for four different pulse durations. 
Simulations with pump pulse durations of 48fs and over never 
reached a steady state, and so are not included here. The first 
point to note in fig. |35]is that the SVEA results are not identi- 
cal in all frames, even though the scaling procedure in section 
nil Al is designed to make them, as far as possible, indepen- 
dent of pulse duration. The reason is that, as noted in section 
nil Al the dispersion scales in a different way to the group time 
delay, and so is not correctly compensated by eqn. ( I20> . 

A second rather surprising feature is that we might expect 
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FIG. 31: Dispersion Scaled OPO: Time domain representation of 
the GFEA modulus-squared of the pulse envelopes, for a range of in- 
jected pump pulse durations: 6-192fs. (bottom to top) pump, signal, 
and idler. 



the GFEA results to tend to the SVEA as pulse length in- 
creases, but this is not evident from the graphs. The expla- 
nation for this is that the steady state of the OPO can change 
suddenly as the parameters are varied. This property is high- 
lighted in fig. 1^ which shows the GFEA signal pulse profile 
for pulse durations from 6fs to 192fs; the sudden adjustment 
of the GFEA when moving from 36fs to 48fs takes it close to 
the SVEA, and the difference between the two gradually dis- 
appears as the pulse duration is increased further (see Fig. |33 
). Note that the scaling procedure used for fig. |^is an ex- 
tension of eqn. (|20} in that the dispersion term is also scaled, 
making the SVEA results completely independent of pulse du- 
ration. 

Spectral profiles corresponding to the temporal profiles of 
fig- and|36]are shown in fig. |^ The spectral shape for 
each field is similar across all pulse durations, with a pulse of 
double the (time) width naturally having half the bandwidth. 




FIG. 32: Dispersion Scaled OPO: Time domain representation of the 
GFEA amplitude phase, for a range of injected pump pulse durations: 
6-192fs. (bottom to top) pump, signal, and idler. 



Notice that the pump and signal spectra in the 6fs frame are 
close to overlapping, which indicates that the separation of the 
total EM field into distinct pump, signal, and idler components 
is becoming a questionable assumption. 

Inclusion of the caiTier wave in the results raises some 
quite subtle issues that need careful consideration. It must be 
stressed again that the carrier drops out of the analysis leading 
to eqn. (^3. The envelope description is therefore complete, 
although the phases of two of the three envelope functions 
can be changed by arbitrary constants without any effect on 
the computations apart from an appropriate adjustment in the 
phase of the third envelope. For instance, if the phases of the 
pump and signal envelopes are changed by A(|)p and A(|)s, the 
phase of the idler envelope is changed by A(|)/ = A(|)p — Ac]), . 
Adjustments of this kind show up in the results only if graphs 
of the complete electric field profiles, including the carrier 
waves, are displayed, as in fig. |38] If the simulations in that 
figure were re-run with differing envelope phases, this would 
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FIG. 33: Dispersion Scaled OPO: Time domain representation of 
the difference between the modulus-squared of the pulse envelopes 
between SVEA and GFEA theories, for a range of injected pump 
pulse durations: 6-192fs. (bottom to top) pump, signal, and idler. 
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FIG. 34: Dispersion Scaled OPO: Maximum difference 5 between 
GFEA and SVEA simulations over the middle quarter of the scaled 
T range, on a logjg scale. 
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FIG. 35: Scaled OPO: Time domain representation of the modulus- 
squared of the pulse envelopes, for a range of injected pump pulse 
durations: 6fs (top), 9, 12, 18, 24, 36fs (bottom). For each sub- 
figure, the curves compare (bottom to top) pump, signal, and idler 
for the SVEA simulations ( ) and GFEA ones ( ). 
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FIG. 36: Scaled OPO: As above in fig. |35l for 6fs (top), 9, 12, 
18, 24, 36fs (bottom); but using the same time window width for all 
pulse widths. 
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FIG. 37: Scaled OPO: Frequency domain representation of the 
modulus-squared of the pulse envelopes, for pump pulse durations 
of 6, 12, 18, 24, and 36fs. For each sub-figure, the curves compare 
(bottom to top) pump, signal, and idler for the SVEA simulations 
( ) and GFEA ones ( ). 
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be reflected in temporal displacements of the carrier-like os- 
cillations beneath the envelopes. 

A further interesting feature is that, while the moduli of the 
pulse enveleopes may have stabilised in a simulation, the en- 
velope phases can (and usually do) change from pass to pass; 
this process continues indefinitely, so a movie made up of 
frames from successive transits would show the pump, signal, 
and idler electric field oscillations drifting across underneath 
the respective steady envelope profiles. 

Note that although it is common to talk of "shifts in carrier 
phase", in the theory the carrier phase is specified by the ini- 
tial conditions, and does not change. The effects usually as- 
cribed to "shifts in carrier phase" instead manifest themselves 
as a shift in the phase of the complex envelope. 

The different models discussed in this paper give signifi- 
cantly different results for the pass-to-pass phase drift. Fig. 
|39lshows the phase change for the signal pulses as a function 
of pulse length for the SVEA, SEWA, and GFEA; note that the 
SEWA and GFEA results are similar to each other, while the 
(less accurate) SVEA exhibits a very different dependence. 

The reference point used in calculating the phase drift is at 
the maximum amplitude of the envelope of the signal pulse, 
which is in fact not necessarily at the point of maximum elec- 
tric field. This is a good choice for our purposes because it 
does not move between passes once a steady state is estab- 
lished. Although these phase drifts are quite small, discrepan- 
cies between the SVEA and GFEA will quickly accumulate. 

Analysis shows that the phase of the envelope at this ref- 
erence changes by a fixed amount from one pass to the next. 
Since this phase drift tends to be small, a slow evolution of 
the E field profile can be seen between passes in any given 
simulation. Further, the accumulated difference between the 
less accurate SVEA prediction and the GFEA can lead to large 
differences in the predicted electric field profiles, even if the 
envelope profiles happened to be similar. Note that the idler 
pulse envelope phase drift is in the opposite direction to the 
signal drift, since the pump pulse arrives with the same enve- 
lope phase at the beginning of each pass. 

The details of fig. |39]are less important than its message. 
Getting an output of either an identical idler pulse from pass to 
pass or at least a well understood pulse shape progression can 
be important, and this will be achieved using accurate models 
or by experiment based on accurate models. We see that even 
for relatively long pulses the SVEA theory predicts a different 
phase drift to that of the SEWA/GFEA theory; and that the 
trend is different even when using a sequence of parameters 
designed to minimise the differences between simulations of 
different pulse lengths. 



C. Synchronisation and Pliase Mismatcli 
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The phase drift varies as a function of the phase mismatch 
and the cavity synchronisation: the data is shown in Figs. |^ 
im The glitches in the data are caused by occasional diffi- 
culties in evaluating the phase of the envelopes and the wide 
variations that can occur for near-zero fields - not all the data 



FIG. 38: Scaled OPO: Time domain representation of the electric 
fields of the pulse, for pump pulse durations of 6, 12, 18, 14, and 

36fs. For each sub-figure, the solid curves ( ) compare (bottom 

to top) pump, signal, and idler for the GFEA simulations, for 6fs and 

9fs the SVEA fields are also indicated ( ). The phases are chosen 

so that the maximum excursion of the signal envelope is purely real 
valued, and the idler phase is chosen so that (j)., + (]),■ = i^p. 
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FIG. 39: Scaled OPO: Pass-to-pass phase drift for a range of param- 
eters, comparing SVEA ( ), GFEA ( ), and SEWA ( ) 

simulation results. The differences are taken between the phase at 
the peak of the modulus-squared of the envelopes at the end of one 
pass of the signal pulse and the next. 



points represent good steady states. 
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D. Fixed Length Crystal 




Here we present results for a more realistic case of a fixed 
length crystal (1000/jm) and a range of pump pulse durations 
(6 - 96fs); with all parameters being kept constant, notably a 
fixed pulse energy. Since the crystal is relatively long, the high 
dispersion comes into play, ensuring that even for rather short 
pump pulses (e.g. 6fs), the signal and idler become relatively 
broad 200fs) - thus the role of few-cycle effects should be 
relatively small. However, inspection of the envelope profiles 
for the SVEA and GFEA models still show small but notice- 
able differences, for both perfect and imperfect phase match- 
ing, even for pump pulses as long as 96fs. 

In fig. 1^ we see how the inter-pass phase difference 
for the signal (and hence idler) for a fixed 1000/jm crystal 
with decreasing pump pulse widths and adjusted time offsets, 
but fixed pulse energy. Difference between the SVEA pre- 
dicted phase drift and the GFEA are noteable for perfect phase 
matching (Afe = 0), but when this is no longer exact, the two 
predictions rapidly become similar, as now both a long crystal 
and imperfect phase matching both act against the generation 
of clear few cycle effects. 
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VI. CONCLUSIONS 

We have presented a new and more complete derivation of 
how the envelopes of extremely short optical pulses evolve in 
nonlinear interactions. We have compared the results of our 
new (GFEA) model to those of the traditional slowly varying 
envelope approximation (SVEA) using a scaling procedure to 
distinguish specific few-cycle effects from other phenomena 
caused by changing pulse duration. It should be noted that the 
SVEA becomes inadequate whenever the envelope changes 
rapidly within a few carrier periods. Strictly speaking, a few- 
cycle pulse is not required, because a steep edge within a 
longer pulse also fulfils the conditions. 

The effect of the extra "few-cycle" terms in the GFEA evo- 




FIG. 40: Scaled OPO: GFEA phase mismatch vs sync 
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FIG. 41: Scaled OPO: SVEA phase mismatch vs sync 
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FIG. 42: 1000/jm OPO: Pass-to-pass phase drifts for a differing 
pulse lengths but a fixed crystal length of lOOOjUm and equal pulse 

energies, comparing SVEA ( ), SEWA ( ), and GFEA ( ) 

simulation results. The differences are taken between the phase at 
the peak of the modulus-squared of the envelopes at the end of one 
pass of the signal pulse and the next. Each graph is for a different 
phase mismatch Ak in the periodic poling of the crystal. 



lution equation is to add a phase distortion to the nonlinear 
polarization term, which then imposes itself on the pulse en- 
velopes. This is demonstrated by our single-pass optical para- 
metric amplifier NPA model where, whilst the SVEA model 
is insensitive to pulse length, the GFEA theory shows clear 
changes as the pulses get shorter and contain fewer optical 
cycles. 

Further, when we studied the highly sensitive de- 
amplification case (i.e. NPD), we saw dramatic differences 
between the SVEA and GFEA simulations even outside the 
few-cycle regime. These arose from the phase distorting ef- 
fects of the few-cycle terms in the theory disrupting the exact 
phase relationships needed for de-amplification. While the 
absolute size of these differences do depend on the chosen pa- 
rameters of crystal length, pulse energy, and so on, they will 
always get dramatically larger for shorter pulses. 

On the other hand, the repetitive cycling nature of the op- 
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tical parametric oscillator (OPO) produces more complicated 
and subtle dynamics; small changes in parameter values can, 
for instance, cause sudden changes in the steady state fields. 
It is therefore no surprise that comparison of the results pre- 
dicted by the different models is less straightforward in the 
OPO case. The new model certainly produces differences in 
the pulse envelopes as well as the phases, although the way in 
which the GFEA tends to the S VEA in the long-pulse limit has 
some interesting features. The two models also predict differ- 
ent results for the pass-to-pass phase drift of OPO pulses, and 
this implies significant differences in the electric field struc- 
tures. In both cases, the carrier wave moves under the enve- 
lope from one transit to the next, but by different amounts. 



'Dr. Paul. Kinsler ©physics . org! 
|http://www.kinsler.org/physicsJ 

It can be useful to regard these "few-cycle" effects as 
adding a phase twist to the envelope evolution. This then 
shows itself most clearly in the pass-to-pass phase drift in the 
signal and idler phases (see Fig. I39> . Finally, although these 
effects are usually described as "few-cycle effects", they are 
more accurately described as "finite pulse length" effects. Al- 
though for the parameter ranges needed to describe typical 
nonlinear crystals, it is only in a few-cycle regime where the 
effects are easily visible in the pulse intensity profiles, these 
twisting effects on the phase structure of the pulse do occur, 
and could be seen in a many cycle regime. 
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